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Abstract 

We consider stochastically quantized self-interacting scalar fields as suitable models to generate 
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chaotic dynamics conjugated to a Bernoulli shift in fictitious time, and the right amount of vacuum 
energy density can be generated without fine tuning. It is numerically observed that the scalar 
field dynamics distinguishes fundamental parameters such as the electroweak and strong coupling 
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possible solutions to the cosmological coincidence problem, as well as to the problem of uniqueness 
of vacua. 

PACS numbers: 98.80.-k, 03.70.+k, 05.45.Jn 


*Permanent address: School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, 
LondonEl 4NS, UK; Electronic address: c.beck@qmul.ac.uk, URL: http://www.maths.qmul.ac.uk/~beck 


1 



I. INTRODUCTION 


There is by now convincing observational evidence that the universe is currently in a 


phase of accelerated expansion |l|, 


The favored explanation for this behavior is the 


existence of vacuum energy or, in a more general setting, of dark energy. The observations 
suggest that the universe currently consist of approximately 73 % dark energy, 23 % dark 
matter, and 4% ordinary matter [3j]. The nature and origin of the dominating dark energy 
component is not understood, and many different models co-exist. The simplest models 
associate dark energy with the vacuum energy of some unknown self-interacting scalar field, 

In quintessence models slowly 
String theory 


Q. 


whose potential energy yields a cosmological constant | 
evolving scalar fields with a nontrivial equation of state are considered j 
also yields possible candidates of scalar fields who might generate dark energy, in form of 
run-awav dilatons and moduli fields |f|. Various exotic forms of matter such as phantom 
matter |tJ and Born-Infeld quantum condensates are currently being discussed. For some 
superstring cosmology ideas related to small cosmological constants, see also jj|. 

When trying to formulate a suitable model for dark energy, at least two unsolved funda¬ 
mental problems arise: 

1. The cosmological constant problem. Why is the observed vacuum energy density 
so small, as compared to typical predictions of particle physics models? From electroweak 
symmetry breaking via the Higgs mechanism one obtains a vacuum energy density prediction 
that is too large by a factor 10 55 as compared to the currently observed value. Spontaneous 
symmetry breaking in GUT models is even worse, it yields a discrepancy by a factor 10 111 . 

2. The cosmological coincidence problem. Why is the order of magnitude of the currently 
observed vacuum energy density the same as that of the matter density? A true cosmological 
constant stays constant during the expansion of the universe, whereas the matter energy 
density decreases with a -3 , where a(t) is the scale factor in the Robertson Walker metric. It 
looks like a very strange coincidence that right now we live at an epoch where the vacuum 
energy density and matter density have the same order of magnitude, if during the evolution 
of the universe one is constant and the other one decreases as a(t)~ 3 . 

To this list one may add yet another fundamental problem, which we may call 

3. The uniqueness problem. String theory allows for an enourmous amount of possible 
vacua after compactification. In each of these states the fundamental constants of nature 


2 


can take on different possible values. But what is the mechanism that selects out of these 
infinitely many possibilities the physically relevant vacuum state, with its associated funda¬ 
mental constants that give rise to a universe of the type we know it (that ultimately even 
enabled the development of life)? Relating the answer purely to an anthropic principle seems 
unsatisfactory. 

In this paper we consider a new model for dark energy which, as compared to other 
models, is rather conservative. It just associates dark energy with self-interacting scalar 
fields corresponding to a </? 4 -theory, which is second quantized. However, the fundamental 
difference to previous approaches is that these fields are very strongly (rather than weakly) 
self-interacting, and that 2nd quantization effects play an important role. We will use as the 
relevant method to quantize the scalar fields the stochastic quantization method introduced 
by Parisi and Wu |10|. In the fictitious time variable of this approach, the fields will turn 
out to perform rapid deterministic chaotic oscillations, due to the fact that we consider not 
a weakly but a very strongly self-interacting field. This chaotic behavior is a new effect not 
present in any classical treatment. It is generally well known that chaos plays an important 


role in general relativity [llj], quantum field theories 1 
The main result of our consideration is that the chaotic field theories considered naturally 
generate a small cosmological constant and have the scope to offer simultaneous solutions 
to the cosmological coincidence and uniqueness problem. 

Our physical interpretation is to associate the chaotic behaviour of the scalar fields with 
tiny vacuum fluctuations which are allowed within the bounds set by the uncertainty relation, 
due to the finite age of the universe. This interpretation naturally leads to the right amount 
of dark energy density being generated, and fine tuning can be avoided. The chaotic fields 
(presently) have a classical equation of state close to w = — 1, and can thus account for the 
accelerated expansion of the universe. However, during the early evolution of the universe 
they behave in a different way: They effectively track radiation and matter. This property 
will help to avoid the cosmological coincidence problem. 

The chaotic model also contains an interesting symmetry between gravitational and gauge 
couplings. In our model the role of a metric for the 5th coordinate (the fictitious time) is 
taken over by dimensionless coupling constants which are given by the ratio of the fictious 
time lattice constant and physical time lattice constant squared (both lattice constants can 
still go to zero, just their ratio is fixed). These coupling constants do not occur in any 


* QQQ, 


and string theories 


Q. 
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classical treatment but are entirely a consequence of our second quantized treatment. The 
vacuum energy generated depends on these couplings in a non-trivial way. The physical 
significance of our model is illustrated by the fact that we numerically observe the vacuum 
energy to have local minima for coupling constants that numerically coincide with running 
electroweak coupling strengths, evaluated at the known fermionic mass scales, as well as 
running strong coupling constants evaluated at the known bosonic mass scales. This numer¬ 
ical observation, previously reported in Q, is now embedded into a cosmological context. 
The role of the chaotic fields in the universe can be understood in the sense that they are 
responsible for fixing and stabilizing fundamental parameters as local minima in the dark 
energy landscape. This is somewhat similar to the role the dilaton field plays in string theory 
after supersymmetry breaking. 

Our numerical discovery of local minima that coincide with known standard model cou¬ 
pling constants makes it very unlikely that there are different universes with different fun¬ 
damental parameters. In fact, the numerical results provide strong evidence that there is a 
unique vacuum state of the universe that possesses minimum vacuum energy precisely for 
the known set of standard model parameters. 

This paper is organized as follows. In section 2 we show how a second-quantized scalar 
field dynamics can degenerate to a chaotic dynamics in fictitious time. Our main example is 
a chaotic y? 4 -theory leading to 3rd order Tchebyscheff maps, which is dealt with in section 
3. In section 4 we present a physical interpretation of the chaotic dynamics using the 
uncertainty relation, which in a natural way fixes the order of magnitude of the vacuum 
energy density to be generated. Section 5 deals with energy, pressure and classical equation 
of state of the chaotic fields. In section 6 we consider the Einstein equations associated with 
our model and discuss a possible way to avoid the cosmological coincidence problem. Section 
7 yields a prediction for the current ratio of matter energy density and dark energy density 
to the critical energy density. In section 8 we describe how local minima of the dark energy 
landscape generated by the chaotic fields can fix the fundamental parameters. Finally, in 
section 9 we discuss spontaneous symmetry breaking phenomena for the chaotic fields. 
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II. STOCHASTIC QUANTIZATION OF STRONGLY SELF-INTERACTING 
SCALAR FIELDS 


Let us consider a self-interacting scalar field <p in Robertson-Walker metric. For a com¬ 
plete theory describing all quantum mechanical fluctuations we need to second-quantize it. 
This can be done via stochastic quantization. In the ParishWu approach of stochastic quan¬ 
tization one considers a stochastic differential equation evolving in a fictitious time variable 
s, the drift term being given by the classical field equation [lOj] . Quantum mechanical ex¬ 
pectations correspond to expectations with respect to the generated stochastic processes in 
the limit s oo. The fictitious time s is different from the physical time t, it is just a 
helpful fifth coordinate to do 2nd quantization. Neglecting spatial gradients the field ip is a 
function of physical time t and fictitious time s. The 2nd quantized equation of motion is 



( 1 ) 


where H is the Hubble parameter, V is the potential under consideration and L(s,t ) is 
Gaussian white noise, ^-correlated both in .s and t. For e.g. a numerical simulation we may 
discretize eq. (HD using 


( 2 ) 

( 3 ) 


s 


nr 


t = iS, 


where n and i are integers and r is a fictitious time lattice constant, <5 is a physical time 
lattice constant. The continuum limit requires r —» 0, 8 —> 0, but we will later argue that 
it makes physical sense to keep small but finite lattice constants of the order of the Planck 
length. We obtain 



This can be written as the following recurrence relation for the field tp l n 


V>Ul = (!-«) | <P l n + j + - <Pn ^ + f (^n +1 + <Pn ^ + T ' noise > ( 5 ) 


where a dimensionless coupling constant a is introduced as 


2 r 


(6) 
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We also introduce a dimensionless field variable <b l n by writing p l n = , where p ma x 

is some (so far) arbitrary energy scale. The above scalar field dynamics is equivalent to a 
spatially extended dynamical system (a coupled map lattice jlf| ) of the form 

U = (1 - a)T(& n ) + | HSa^ - ^ _1 ) + |« +1 + ^ _1 ) + r • noise, (7) 

where the local map T is given by 

r($) = $ + a-/ M ' 

Here the prime means 

,_d__ _L _ d_ 

d<p p max <9$ ’ 

Note that a symmetric diffusively coupled map lattice of the form 

& n+1 = (1 - a)T(& n ) + -($** + SjT 1 ) + T ■ noise ( 10 ) 

is obtained if H5 « l, equivalent to 

5«H~\ (11) 

meaning that the physical time lattice constant S is much smaller than the age of the 

universe. In this case the term proportional to H in eq. m can be neglected. The local map 
T depends on the potential under consideration. Since we restrict ourselves to real scalar 
fields <p,T is a 1-dimensional map. 

The main result of our consideration is that iteration of a coupled map lattice of the 
form m with a given map T has physical meaning: It means that one is considering the 
second-quantized dynamics of a self-interacting real scalar field p with a force V' given by 

— ( 12 ) 

Integration yields 

V(<p) — 1 a j- 7^ 2 +Pmax J dp T + const. (13) 

In terms of the dimensionless field $ this can be written as 

V(p) = ^—-^p 2 max j-l$ 2 + J d$T($)| + const. (14) 


( 8 ) 

(9) 
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III. CHAOTIC (/-THEORY 


An interesting observation is the following one. The lattice constant r of fictitious time 
should be small, in order to approximate the continuum theory, which is ordinary quantum 
held theory. If r is small one naively expects the map T given by eq. (JBJ) to be close to the 
identity for finite forces V', since rV'/p max is small. What about, however, very strong forces 
V' due to very strongly self-interacting fields? If p max /r is of the same order of magnitude 
as V' then a nontrivial map T can arise. In particular, this map may even exhibit chaotic 
behaviour. 

As a distinguished example of a /-theory generating strongest possible chaotic behaviour, 
let us consider the map 

$ n+1 = T_ 3 ($„) = -4$^ + 3$ n (15) 

on the interval $ 6 [—1,1]. T_ 3 is the negative third-order Tchebyscheff map, a standard 
example of a map exhibiting strongly chaotic behaviour. It is conjugated to a Bernoulli 
shift, thus generating the strongest possible stochastic behavior possible for a smooth low¬ 
dimensional deterministic dynamical system. The corresponding potential is given by 

V_ 3 (/ = —jV - Z| # const , (16) 

or, in terms of the dimensionless field <f>, 

V_ 3 (Z = -— Pmaxi® 2 ~ ^ 4 ) + Const- (17) 

T 

Apparently, starting from this potential we obtain by second quantization a field ip that 
rapidly fluctuates in fictitious time on some finite interval, provided that initially ip 0 e 
[—Pmax 1 Pmax]- The small noise term in eq. (fTUl) can be neglected as compared to the deter¬ 
ministic chaotic fluctuations of the field. 

Of physical relevance are the expectations of suitable observables with respect to the 
ergodic chaotic dynamics. For example, the expectation (VL. 3 (/) of the potential is a 
possible candidate for vacuum energy in our universe. One obtains 

(V_ 3 {/) = -—- P 2 max (($ 2 ) - ($ 4 )) + con st- (18) 

T 

For uncoupled Tchebyscheff maps (a = 0), expectations of any observable A can be evaluated 
as the ergodic average 

(*) = f l (is) 
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with the natural invariant measure being given by 


d<3> 

d ( 20 ) 

7TVl - $ 2 

(see any textbook on chaotic dynamics, e.g. Q). This measure describes the probability 
distribution of the iterates under long-term iteration. From eq. (12Uj) one obtains (<f> 2 ) = \ 
and (<h 4 ) = |, thus 

(V-M) = + const . ( 21 ) 

Alternatively, we may consider the positive Tchebyscheff map T 3 (<&) = 4$ 3 — 3$. This 
basically exhibits the same dynamics as T_ 3 , up to a sign. Repeating the same calculation 
we obtain 

V 3 [<fl) = j -2 ^ 2 + 1 P 4 j + const (22) 

and 

V 3 (ip) = ^—^P 2 max (- 2$ 2 + $ 4 )- (23) 

T 

For the expectation of the vacuum energy one gets 

(V^p)) = -- -Pmax(- 2(® 2 ) + (® 4 )) + const, (24) 

r 

which for a = 0 reduces to 

{VM) = - + cans t. (25) 

Symmetry considerations between T_ 3 and T 3 suggest to take the additive constant const 
as 

(mist. - -j~-—— PmaxTii® 2 )- (26) 

r 2 

In this case one obtains the fully symmetric equation 

(V a (v)) = {-|<® 2 } + (4> 4 )} , (27) 

which for oi —> 0 reduces to 

(V ±3 ( V )) = (_?) . (28) 

The simplest model for dark energy in the universe, as generated by a chaotic <p 4 -theory, 
would be to identify fp^ax/ 7 " = Pa, the constant vacuum energy density corresponding 



to a classical cosmological constant A, which stays constant during the expansion of the 
universe. This is certainly a possible simple model. On the other hand, such an approach 
would neither solve the cosmological constant nor the cosmological coincidence problem. 
For this reason we will turn to a more advanced model in the following sections, which will 
naturally produces the right amount of vacuum energy density in the universe. 

Before proceeding to this model, let us provide some general comments on the physical 
meaning of the parameter r. The vacuum energy generated by the chaotic fields is inversely 
proportional to r (see eq. Strict equivalence of the stochastic quantization method 

with quantum field theory requires the continuum limit r —► 0. In this limit eq. (J2%|) can still 
generate a finite amount of vacuum energy, provided both r —> 0 and p max —> 0 such that 
the ratio p^ ox /r stays finite. In fact, many of the results presented in this paper depend 
only on the ratio p 2 iax /T and not on the individual values attributed to r and p max ■ 

From the viewpoint of ordinary quantum held theory, the limit r —► 0 and p ma x 0 
means that one formally considers a potential V(<p) — p 2 p 2 + Ap 4 for which both potential 
parameters p 2 ~ r _1 and A ~ r _1 p m 2 x diverge (see eq. dm moreover the held ip lives on an 
inhnitely small support [—p m ax, Pmax]- Clearly, this is a very singular type of quantum held 
theory, which in principle can also be studied by other means than stochastic quantization, 
though perturbation theory will not be applicable. The advantage of our formulation in 
terms of stochastic quantization is that for the dimensionless chaotically evolving held $ the 
potential remains hnite (see, e.g., eq. (|28|) h If there is no hctitious time, then the parameter 
r enters in form of the (singular) potential parameters p 2 ~ r _1 and A ~ r _1 p m ax- I 11 the 
next section we will argue that on physical grounds it makes sense to consider very small 
but hnite r. 

IV. REPRODUCING THE CURRENTLY MEASURED DARK ENERGY DEN¬ 
SITY 

To obtain quantitative statements on the dark energy density as generated by some 
chaotically evolving held <p, let us hx the free parameters r and p max by some physical 
arguments. Let us start with the parameter r. It is the lattice constant of hctitious time s 
and has dimension GeV~ 2 . Ordinary stochastic quantization based on Gaussian white noise 
requires the continuum limit r —> 0. But since quantum held theory runs into difficulties 
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at the Planck scale mpi and is expected to be replaced by a more advanced theory at this 
scale, it is most reasonable to take the small but finite value 

r ~ m^. (29) 


Next, consider the parameter p max • It has dimension GeV and describes the maximum 
energy scale of our rapidly fluctuating scalar fields ip, who take on values on the finite interval 
[—Pmax,Pmax)- A natural value of p max follows if one associates the rapidly fluctuating chaotic 
fields p % n with vacuum fluctuations that are allowed due to the uncertainty relation 

AEAt = Q(h). (30) 


Taking At ~ t of the order of the age of the universe, a corresponding energy uncertainty 
A E arises. This A E is very large for a very young universe, and then decreases to extremely 
small values for the current age of the universe of about 13.7 Gyr. Any finite age t of the 
universe implies that spontaneous vacuum fluctuations with energies of order A E ~ f -1 can 
occur. It is physically plausible to identify these energy fluctuations A E with the rapidly 
fluctuating chaotic fields ip — p m ax®h, since both A E and ip live on a finite interval, and 
both fluctuate in an unpredictable way. The uncertainty relation (HU1) together with At ~ t 
implies (in units where h — c = 1) 

(31) 

In this way the energy scale p max occurring in the chaotic field theories is most naturally 
identified with the inverse age of the universe. However, note that this quantum mechanical 
interpretation in terms of vacuum fluctuations requires a chaotic map T, since some regularly 
evolving cannot be associated with fluctuations at all. 

It is remarkable that by taking eq. (EH and eq. (ED together, the right amount of vacuum 
energy follows without any fine tuning. One has for generic chaotic maps T 

(V(*>)> ~ ~ H 2 m 2 pl . (32) 


since t 1 ~ H. Moreover, 


ff 2 = —ft, 

3 P 


where p c denotes the critical density of a flat universe and G 
constant. Combining eq. (T721) and (I33D one obtains 


(33) 

- mp 2 is the gravitational 


(V(v))~Pc, 


(34) 
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as required and confirmed by current astronomical observations. Our simple physical inter¬ 
pretation, namely to interpret the chaotic fluctuations as vacuum fluctuations allowed due 
to the finite age of the universe, thus yields the right order of magnitude of dark energy 
density. 

In general, eq. m only yields order-of magnitude estimates. Nevertheless, it is in¬ 
structive to work out some concrete numbers, based on simple model assumptions. For 
example, we may assume that the entire vacuum energy of the universe is due to one 
chaotic held described by V_ 3 . The current age of the universe is t 0 — (13.7 ± 0.2) Gyr 
= (4.32 ± 0.06) • 10 17 s p. Using an uncertainty relation of the form AEAt = h/2 we 
get p max = 1/(21 0 ) = (7.62 ± 0.08) • 10~ 43 GeV. Choosing r = nrripf, where k is some 
dimensionless number of O(l), we get 

{V(<p)) = = (3.19 ± 0.05) • KT 47 ^ 1 GeV 4 (35) 

The current observational estimate of dark energy density in the universe is jii| 

p° bs = (2.9 ± 0.3) • 10“ 47 0eU 4 , (36) 

which is consistent with k ~ 1. If the observed dark energy in the universe is produced by 
our chaotic theory, then the measured data imply 

k= 1.10 ±0.10. (37) 

Once again, this estimate is based on concrete model assumption. In general, we do do not 
know the precise values of the proportionality constants in eq. (I2U1) and (ETT)) . moreover we 
do not know how many different chaotic fields may contribute to the dark energy of the 
universe. 

Of course, the actual properties of the chaotic dark energy component depend on the 
classical equation of state of the chaotic fields, which will be dealt with in the next section. 

V. ENERGY DENSITY, PRESSURE, AND EQUATION OF STATE 

The kinetic energy term of our chaotic fields is given by 

Tkin =2 (di^) ■ (38) 
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Discretized with lattice constant <5 we obtain for the expectation of T kin 


(T tm ) = ^((41-4-r 1 ) 2 ) 

= ^ia«4> 2 > - (KK')) (39) 

In particular, for a —» 0 the expectation of kinetic energy vanishes, and a universe mainly 
filled with such a field is vacuum energy dominated. 

In general, the expectation of the total energy density {p) is given by 

(p) - (T kin ) + (V) (40) 

and the expectation of the pressure by 

(p) - (T ki n) - (V }. (41) 


For the map T_ 3 one obtains 

M = ^ {f ((4> 2 > - (4i,4r')) + (1 - a)(|<4> 2 > - <4» 4 ))} 

<P> = ^ {f «4> 2 > - («“'» - (1 - <*)(§<4> 2 } - (* 4 ))} 

where the additive constant of the self-interacting potential is fixed by the 
sideration of section 3. The above equations yield the equation of state 


(42) 

, (43) 

symmetry con- 

(44) 


which varies as a function of the coupling a in a nontrivial way. 

For a = 0, the equation of state of our fields is w = — 1, since the expectation of kinetic 
energy vanishes (all expectations should be interpreted as quantum mechanical expectations 
with respect to second quantization). For small a, w is close to —1. It should be clear that 
although our fields fluctuate rapidly in both physical and fictitious time, these fluctuations 
are averaged away when doing the quantum mechanical expectations. Thus the classical 
picture that arises out of this 2nd quantized rapidly fluctuating model is a very homogeneous 
field. 

The expectations in eqs. (|32j) . (pEIj) are easily numerically calculated by long-term iterating 
the coupled map lattice for random initial conditions and averaging over all i and n. We 
used lattices of size 10000 with periodic boundary conditions. The result for the equation 
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of state w(a) is displayed in Fig. 1. For small a, w grows approximately in a linear way. 
It monotonously increases from w = —1 for a = 0 to w = +1 for a = 1, up to a wiggle 
at a « 0.12. Fig. 2 shows the corresponding (classical) energy density and pressure of the 
field. To account for the currently observed dark energy in the universe, most chaotic fields 
must have a coupling a smaller than about 0.04. Larger a are ruled out by the observations 
providing evidence for w < —0.78 jij. 

VI. EINSTEIN EQUATIONS AND DYNAMICAL EVOLUTION 


Let us now consider the Einstein equations for a homogeneous and isotropic universe that 
consists of three different components: matter, radiation, and chaotic fields. These three 
components are labeled by the indices m, r, </?, respectively. One has 

(45) 

(46) 


H 2 — y- J — -7rG(/%, + p m + Pr) 
a = ~^ 7rG (' C V + 3 2V + Pm + Pr + 3 Pr)- 


Here pj denotes the (classical) energy density of component j, and p 3 the pressure. For sim¬ 
pler notation we omit the expectation values (•••). The equation of state of each component 
is Wj = Pj/Pj- As it is well known, one has for matter w m = 0, for radiation w r = 1/3, 
whereas the equation of state w v of the chaotic fields ip depends on the coupling a (see 
Fig. 1). 

The Einstein (or Friedmann) equations are usually supplemented by the assumption of 
conservation of energy for each species j, 


Pj = “3 H( Pj +pj). 


(47) 


These equations can be derived from the Einstein equations under the assumption of adia¬ 
batic expansion, i.e. one assumes that no entropy is produced. 

For a universe dominated by a species j with constant equation of state Wj = Pj/pj 
eq. (jUD leads to 

Pj ~ <r 3(1+ "V (48) 

We obtain the well-known result that for matter p m ~ a -3 , for radiation p r ~ a -4 , whereas 
for true classical vacuum energy (a cosmological constant A) with w = — 1 one has no 
dependence on a at all, p\ — const. 


13 


For the chaotic fields, energy conservation is a non-trivial issue, for the following reasons: 
i) These fields model vacuum fluctuations, and vacuum fluctuations by definition do not 
conserve energy ii) The chaotic field dynamics (as any chaotic dynamics) constantly produces 
entropy in fictitious time, whereas the Friedmann equations (containing no fictitious time) 
describe an adiabatic expansion of the universe iii) The coupling constant a and hence the 
equation of state of the chaotic field may change in time iv) There may be an entire spectrum 
of chaotic fields with different couplings a who can interact with each other v) The chaotic 
fields may interact with dark matter. 

In the following, we want to restrict ourselves to a simple scenario where the coupling a 
is small and where there is energy conservation of quantum mechanical expectations in full 
agreement with the Friedmann equations. Two interesting possibilities arise in this context: 

a) a symmetric phase of the universe, where both the negative and positive 
Tchebyscheff dynamics are active. In this case the vacuum energies p~ := (V_ 3 (</?)) and 

:= (b+ 3 (</?)) cancel to zero: 

Pv = Pp + Pip = 0 ( 49 ) 

(a consequence of eq. (1271) 1. Since the total vacuum energy p lf is zero, there is no statement 
from the Einstein (or Friedmann) equations on the time evolution of the absolute value |p±|. 
We are free to postulate the validity of eq. (1221) 

\pt\ = \’y±M)\~^f L ~ r ^ L (so) 

for arbitrary couplings a. 

b) an asymmetric phase of the universe, where only the negative Tchebyscheff dynamics 
T_ 3 is active (or where it dominates as compared to T + s). In this case one has 

P V > 0- (51) 

If there is a chaotic field with coupling constant a that does not interact with other fields 
then energy conservation implies 

p v ~ a~ 3 ^ 1+Wip \ (52) 

where w v — vj :p (a) is the equation of state as a displayed in Fig. 1. Note that chaotic fields 
with rw ¥ ,Ri0(Q!Ri0.25) behave similar as dark matter, whereas chaotic fields with w v « — 1 
(a & 0) act like a cosmological constant (see also Q). 
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We now show that eq. CTl . derived from the uncertainty relation, naturally leads to 
tracking behavior of vacuum energy density. Suppose the universe surrounding the chaotic 
fields is dominated by a species with equation of state w > — 1 (typically matter or radiation), 
then 

a(t) ~ tTO <*=► r 2 ~ a- 3{1+w) . (53) 

Putting this into eq. m we obtain 

l4M043(‘ri)l~«- 3<1+ ”’ ) , (64) 

i.e. the vacuum energy density associated with the chaotic fields decays in the same way 
with a as the density of the dominating species. 

Eq. (JH3J) can help to naturally avoid the cosmological coincidence problem. Consider 
e.g. the following scenario. Initially (say, shortly after inflation) we may have a symmetric 
phase where p r ~ \p±\ » p m . Here p~ denotes the vacuum energy density of the negative 
Tchebyscheff map, which is cancelled by the vacuum energy density = — p~ of the 
positive Tchebyscheff map. Then \p±\ first decays approximately as a -4 , since the universe 
is radiation dominated. At some stage we arrive at p r ~ \p±\ ~ p m , and from then on 
matter dominates over radiation, so that from then on \p±\ decays approximately as a -3 . 
During the late-time evolution of the universe, \p±\ will always stay of the same order of 
magnitude as p m , since both \p±\ and p m decay as a~ 3 . In spite of this, for small enough 
couplings a the chaotic fields have a classical equation of state close to w v — —1, and 
can thus produce the accelerated expansion of the universe via eq. (flfij) . provided there is 
symmetry breaking between T+ 3 and T_ 3 at some late stage of the evolution of the universe. 
A concrete mechanism for this will be discussed in section 9. 

What is our physical interpretation of the chaotic fields in the universe? For a = 0, 
it can be rigorously proved that rescaled deterministic chaotic Tchebyscheff maps can be 
used to generate spatio-temporal Gaussian white noise on a larger scale QQ . In other 
words, on fictitious time scales r' » r and physical time scales S' » 5 the chaotic noise 
just looks like ordinary Gaussian white noise. We may thus couple the chaotic fields ip to 
ordinary standard model fields in order to second quantize the standard model fields, i.e. 
use the chaotic fields as a source of quantization noise. This is the basic idea of the so-called 
’chaotic quantization’ approach Q . The chaotic fields are well embedded in this way and 
since they are just playing the role of quantization noise, we do not expect them to have 
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any disturbing influence on, say, baryogenesis and similar processes in the early universe. In 
this interpretation vacuum energy just arises out of the expectation of a classical potential 
that generates quantization noise. 


VII. PREDICTION OF % AND Q m 


Our approach allows for the prediction of the order of magnitude of cosmological pa¬ 
rameters such as Qp = pylp c and Q m = p m /p c at the present time. Let us start from the 
uncertainty relation in the form 


AEAt = -, 
2 ’ 


(55) 


which implies 


Choosing the time scale r = Km P f we get for a ps 0 

p; = <r_ 3 M> = = 


3 m 2 Pl 1 


(56) 


(57) 


During the radiation dominated period of the universe one has for the energy density of 
radiation 

p r = ^ N(T)T\ (58) 

where T is the temperature and N(T ) is the number of relativistic particle degrees of freedom. 
There is also a relation between time and temperature, namely 


t = 


90 mpi 
32tt 3 N(T) T 2 ' 


(59) 


Putting eq. into © one obtains 


p; = ^ 3 l N ^ T ‘ = l*- < 6 °) 

This equation once again shows that it is reasonable to assume that p r and p~ have the 
same order of magnitude. Since p~ decays in the same way as p r , eq. (JHnj) is valid during 
the entire radiation dominated epoch. Finally p r falls below p m and from then on we have 


7T 

P v ~ ~Pm- 


(61) 
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After symmetry breaking we have p v = « const. This implies a prediction for := 

Pv!Pc at the present time, namely 


% 


_ P<p _ 

Pip + Pm + Pr 


*•/« 

1 +7t/k’ 


(62) 


neglecting p r at the present epoch and assuming that the symmetry breaking took place 
rather recently. In section 4 we saw that the currently observed dark energy density is best 
fitted by the value k= 1.10 ± 0.10. Eq. (1521) yields with this value the prediction 


» 0.74 Q m w 0.26, (63) 

which is consistent with observationsjjj. Again, due to the reasons that were already men¬ 
tioned at the end of section 4, eq. (63) should only be regarded as an order-of-magnitude 
estimate. 


VIII. FIXING FUNDAMENTAL PARAMETERS 

We have seen that chaotic fields can generate the right amount of vacuum energy and 
have the scope to avoid the cosmological constant and coincidence problem. We now show 
that they also offer solutions to the problem of uniqueness of vacua. 

First, let us slightly generalize the chaotic field dynamics (jl Of) to 

4>Ih = (1 - + r 6 (*j l +1 )) (64) 

(we neglect the small noise term). The case a = +1 is called ’diffusive coupling’, the case 
cr = — 1 ’anti-diffusive coupling’. Chaotic fields with 6=1 are called to be of ’type A’ ( 
T x (3>) =: T(4>)), chaotic fields with 6 = 0 to be of ’type B’ (T°(<f>) =: <f>). In Q the chaotic 
fields were called ’chaotic strings’, but this is only a different name for the same dynamics. 
Our derivation in section 2 lead to chaotic fields of B-type with diffusive coupling, but from 
a dynamical systems point of view all 4 degrees of freedom (6 = 0,1, cr = ±1) exist and are 
of relevance. As shown in detail in fl^ , there are two different types of vacuum energies for 
the chaotic fields, namely 
1 . the self energy 

V(a) := ^ (|( 0 2 ) - <# 4 )) ( 65 ) 

and 
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2. the interaction energy 


W(a) := (66) 

2 r 

Eq. (}65j) actually represents the self energy of the map T_ 3 , the self energy of the map T + 3 
has opposite sign and cancels the self energy of T_ 3 in the symmetric phase. 

Basically, the self energy is the expectation of the potential that generates the chaotic 
dynamics in fictitious time, and the interaction energy is the expectation of the potential 
that generates the diffusive coupling in physical time. One may also define a total vacuum 
energy as i? ± (o;) := V(a) ± aW(a), where the — sign corresponds to diffusive and the + 
sign to anti-diffusive coupling. All additive constants are fixed by the postulate of invariance 
of the theory under global and local ^-transformations fb| . In other words, we allow for 
the existence of a symmetric phase. For small a the interaction energy can be neglected as 
compared to the self energy, moreover, the type-A and type-B forms are observed to have 
the same self energy in this limit. 

The central hypothesis of this paper is a symmetry between standard model coupling 
constants and the chaotic field couplings a. We assume that for any dimensionless coupling 
constant a that appears in the standard model of electroweak and strong interactions, there 
is a corresponding chaotic field that is just coupled with this a. The universe then tries to 
reach a state of minimum vacuum energy by adjusting its free parameters in such a way 
that the chaotic fields reach a state of minimum vacuum energy. 

While at first sight this may look like a purely theoretical concept, there is numerical 
evidence that this principle is indeed physically realized. As an example, Fig. 3 shows the 
self energy V(a) = (V -$(</?)) of our chaotic fields of type A with diffusive coupling in the 
low-coupling region. We observe that V(a) has local minima at 


Oi 

= 0.000246(2) 

(67) 

a 2 

= 0.00102(1) 

(68) 

03 

= 0.00220(1) 

(69) 


(ai and a 3 are actually small local minima on top of the hill). 

On the other hand, in the standard model of electroweak interactions the weak coupling 
constant is given by 


®-vieak — ^el " 


(F 3 — Q sin 2 6w) 2 
sin 2 9w cos 2 Ow 


(70) 
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Here Q is the electric charge of the particle (Q = —1 for electrons, Q = 2/3 for w-like 
quarks, Q = —1/3 for d-like quarks), and T 3 is the third component of the isospin (T 3 = 0 
for right-handed particles, T 3 = — \ for e L and d L , T 3 = +~ for u L and u L ). Consider right- 
handed fermions Jr. With sin 1 2 9w = sf — 0.2318 (as experimentally measured Q) and the 
running electric coupling a e i(E ) taken at energy scale E = 3m f we obtain from eq. (fTUj) the 
numerical values 


: ( 3 m d ) 

= 0.000246 

(71) 

: (3m c ) 

= 0.001013 

(72) 

,(3m e ) 

= 0.00220. 

(73) 


There is an amazing numerical coincidence between the local minima a %, c* 2 , « 3 of V (a) and 
the weak coupling constants of f r = Ur, cr , Br, respectively. 

Now regard the fine structure constant a e i and the Weinberg angle sin 2 9 W as a priori 
free parameters. Suppose these parameters would change to slightly different values. Then 
immediately this would produce larger vacuum energy V(a), since we get out of the local 
minima. The system is expected to be driven back to the local minima, and the fundamental 
parameters are stabilized in this way, provided the universe is in an asymmetric phase. 

The above example is only one example of a large number of numerical coincidences ob¬ 
served. In QQ an extensive numerical investigation of self energies, interaction energies, 
and total vacuum energies was performed for the above chaotic field theories. A large num¬ 
ber of amazing numerical coincidences was found Q. These results are described in detail 
in [3, Q, we here only summarize the main results. 


1. The smallest (stable) zeros of the interaction energy W(a) coincide with ru nning 
electroweak coupling constants, evaluated at energies given by the smallest fermionic 
mass scales. Type (A) describes d-quarks and electrons interacting electrically, type 
(B) ?j-quarks and neutrinos interacting weakly. 

2. Local minima of the self energy V(a) coincide with running weak coupling constants 
of right-handed fermions, evaluated at the lightest fermionic mass scales. 

3. Local minima of the total vacuum energy H + (a ) occur at running strong coupling 
constants evaluated at the lightest baryonic energy scales. 
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4. Local minima of the total vacuum energy H (a) occur at running strong couplings 
evaluated at the lightest mesonic energy scales. 


In jl3, |14| also chaotic fields corresponding to 2nd order Tchebyscheff maps were inves¬ 
tigated (see Appendix), and the following numerical coincidences were found: 


1. The smallest (stable) zeros of W{a) coincide with running strong coupling constants 
evaluated at the smallest bosonic mass scales. Type (A) describes the W boson, type 
(B) the Higgs boson. 


2. Local minima of the self energy V ( a) coincide with Yukawa and gravitational couplings 
evaluated at the fermionic mass scales. 

For more details, see HQ- 

All these numerically observed coincidences are not explainable as a random coincidence. 
Rather, they suggest to interpret the coupling constant a of our second-quantized chaotic 
fields ip as a running gauge coupling. We are free to identify a = 2 t/S 2 with a gauge 
coupling, since the occurence of a ratio of lattice constants r and 5 2 is a new effect in our 
2nd quantized discretized theory, and there is no theory of this dimensionless number so 
far, which represents a kind of metric for the 5th coordinate (the fictitious time). So we 
are indeed free to make the hypothesis that a coincides with a running gauge coupling. By 
doing so, we implicitly construct a symmetry between gauge couplings and gravitational 
couplings, since usually the strength of the kinetic term in the action of a field is determined 
by the metric, i.e. gravitational effects, whereas here it is fixed by standard model coupling 
strengths. The chaotic fields appear to select out of the infinitely many vacua allowed by 
string theory the unique ground state that corresponds to the known coupling constants of 
the universe. All free parameters are fixed in the sense that if the fundamental parameters 
(masses, coupling constants, and mixing angles) had different values, larger vacuum energy 
would arise. 


IX. SPONTANEOUS SYMMETRY BREAKING AND CANCELLATION OF UN¬ 
WANTED VACUUM ENERGY 

Chaotic scalar fields not only allow for a simple mechanism to produce dark energy, they 
also yield a simple mechanism to cancel unwanted dark energy. If we assume that both 
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the positive and negative Tchebyscheff dynamics are physically realized, the corresponding 
vacuum energies precisely cancel for symmetry reasons (see eq. (J2HJ))- This symmetry is a 
Z 2 symmetry which is not there for ordinary smoothly evolving scalar fields (where opposite 
potentials lead to unstable or ill-defined behavior). On the other hand, if only the negative 
Tchebyscheff field dynamics is active, or if it dominates, then positive dark energy arises. 
This positive dark energy can drive inflation, fix standard model parameters as local minima 
in the dark energy landscape, and generate late-time acceleration. It is therefore desirable 
to construct a theory that allows for Z 2 symmetry breaking between positive and negative 
Tchebyscheff maps. 

It is clear that in order to fix fundamental parameters with the methods described in the 
previous section, we must have a broken Z 2 symmetry at some stage of the evolution of the 
universe. Indeed, the minimum requirement we need is at least one very early stage of broken 
symmetry, in order to first-time fix the fundamental parameters to the values which make 
the universe work, and another late-time asymptotic state of broken symmetry, in order to 
stabilize the parameters to their known values so that they cannot drift away to other values. 
It is natural to identify the first phase of broken symmetry with the inflationary phase j^j], 
and the other phase of broken symmetry with the late-time state of the universe. Inbetween, 
we may allow for a symmetric state, which has the advantage that nucleosynthesis is not 
spoilt [27]]. 

As a concrete simple model, consider a scalar field a which takes on the value a = 0 
in the symmetric phase and the values a = ±1 in the phase where the Z 2 -symmetry is 
spontaneously broken. The total potential describing the chaotic field dynamics is given by 

V(a, J>_, *+) = + i±^l4s(4> + ), (74) 

where 

V- 8 (*-) = + i(*p) (75) 

is the potential generating the negative Tchebyscheff field dynamics and 

r + 3(4>+) = A=(-24> 2 + + *4 + i(«2)) (76) 

the one generating the positive Tchebyscheff field dynamics in fictitious time (aw 0). In 
the symmetric phase (a = 0) we obtain from eq. (|7ij) 

(V(a, $_,$+)) = 0, (77) 
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wheras in a broken phase with a = — 1 we obtain 

<V>, 4>_, *+)) = (V-s(*-)> = ^(?<4> 2 ) - <4> 4 }) > 0, (78) 

where we have re-labeled $_ = $. 

We assume that the symmetry is first spontaneously broken to a = —1 at the onset 
of inflation. A large amount of positive vacuum energy is generated via eq. (J75|) . since at 
this stage the universe is very young and p max ~ f _1 ~ H. The chaotic fields can help to 
drive inflation, and fundamental parameters are pre-fixed as local minima in the dark energy 
landscape. 

Then, there is a symmetric phase with a = 0. The consideration of section 7 applies 
but the dark energy is suppressed due to symmetry reasons. Big bang nucleosynthesis and 
galaxy formation can go ahead without any problems. Note that during the symmetric 
epoch the fundamental parameters are no longer stabilized as local minima in the dark 
energy landscape. They can drift to slightly different values. This is consistent with the 
experimental findings of a varying fine structure constant Q. 

Finally, there is late-time symmetry breaking to a = —1. This phase is necessary because 
otherwise the fundamental parameters would keep on drifting to different values. By the 
late-time symmetry breaking, the parameters are finally forced back and stabilized at their 
equilibrium values, which were already pre-fixed during inflation. This gives physical sense 
to the role of late-time dark energy in the universe. 

Late-time symmetry breaking could be physically understood as follows. Suppose the 
negative Tschebyscheff dynamics (which generates positive vacuum energy) is always uni¬ 
formly distributed in space. It is a property of empty space-time. At a very early stage, the 
corresponding vacuum energy was strong and may have driven inflation. Then matter and 
radiation is created. Assume that the matter and radiation particles are second-quantized 
by chaotic noise generated by the positive Tchebyscheff map (which generates negative vac¬ 
uum energy). Then, as soon as sufficiently many matter and radiation have been created 
and quantized, the two vacuum energies (h+ 3 ) and (V_ 3 ) can precisely (or almost) cancel, as 
long as matter and radiation are uniformly distributed in space. Inflation may stop in this 
way and we obtain a symmetric phase of the universe after inflation. But at a late stage of 
the evolution of the universe matter clumps into galaxies. Since the negative vacuum energy 
is generated by chaotic quantization noise for each particle it follows the spatial distribution 
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of matter. The positive vacuum energy remains uniformly distributed. Hence after struc¬ 
ture formation there is no spatially uniform cancellation of vacuum energy anymore. Empty 
space has an excess of positive vacuum energy, galaxies are spatial regions with an excess 
of negative vacuum energy due to quantization noise (the negative vacuum energy in the 
galaxy can be partially compensated by positive kinetic terms that arise out of spatial inho- 
mogenities in the galaxy). In this physical interpretation the late-time symmetry breaking 
is related to structure formation. 

X. CONCLUSION 

We have presented a new model for dark energy in the universe. This model is based 
on a rather conservative approach, the assumption of the existence of second quantized 
self-interacting scalar fields described by a <p 4 -theory. However, the main difference is that 
these fields are strongly self-interacting, rather than weakly. When doing 2nd quantization 
using the ParishWu approach, rapidly fluctuating chaotic fields arise. The expectation of 
the underlying potentials yields the currently observed dark energy density. 

The advantage of this new chaotic model is that many of the questions raised in the 
introduction seem to have natural solutions. The cosmological constant problem is avoided, 
in our model the right order of magnitude of vacuum energy is naturally produced if we 
interpret the chaotic dynamics in terms of vacuum fluctuations allowed by the uncertainty 
relation, for a given finite age of the universe. The cosmological coincidence problem is also 
avoided, since in our model the generated dark energy is not constant anymore, but thins out 
with the expansion of the universe in the same way as the energy density of the dominating 
species (matter or radiation). In spite of that, the (classical) equation of state of the chaotic 
component is close to w = —1, and can account for the accelerated expansion of the universe, 
provided there is late-time symmetry breaking. The chaotic fields are physically interpreted 
in terms of vacuum fluctuations. As such they can temporarily violate energy conservation, 
but quantum mechanical expectations are fully compatible with the Friedmann equations. 

The physical relevance of our model is emphasized by the observation of a large number 
of numerical coincidences between local minima in the dark energy landscape and running 
standard model coupling constants evaluated at the known fermionic and bosonic mass 
scales. It thus appears that chaotic fields have the potential to fix and stabilize fundamen- 
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tal parameters and to select the physically relevant vacuum state out of infinitely many 
possibilities. 

Appendix A: General Tchebyscheff maps 

Our approach can be easily generalized to Tchebyscheff maps of arbitrary order N. One 
has Ti($) = <f>, T 2 (<f>) = 2$ 2 — 1, T 3 (<f>) = 4$ 3 — 3$, generally Ty(<3>) = cos(Aarccos$) with 
$ G [— 1,1], A Tchebyscheff map of order N is conjugated to a Bernoulli shift of N symbols, 
it is ergodic and mixing for N > 2. It exhibits the strongest possible chaotic behaviour 
that is possible for a 1-d smooth map, characterized by a minimum sceleton of higher-order 
correlations Q. 

It is useful to consider both positive and negative Tchebyscheff maps and to define 

T-n(Q) := —Tjv( < 3>). (79) 

The behaviour of T_ N under iteration is identical to that of T N up to a sign, the trajectory 
of TLjv differs by a constant sign (N even) or an alternating sign (N odd) from that of T N . 

Eq. (fHl) implies that the maps T N correspond to potentials V N given by 

Vjv(^) = | ~<P 2 +Pmax J dipT N |4 const (80) 

In particular, one obtains for N = ±1, ±2, ±3 

Vtx(<p) - ± + const (81) 

V&fa) = (^-^ 2 ± - Pmax<p)^J + const (82) 

( _ ^ 2 - ± (-^ 2 + + const. (83) 

Of physical relevance are the expectations of these potentials, formed with respect to the 
ergodic dynamics. Since negative and positive Tchebyscheff maps generate essentially the 
same dynamics, up to a sign, any physically relevant expecatation should also be the same 
for T n and T _ N , up to a possible sign. For all N, this symmetry condition fixes the additive 
constant to be 

l Oi l.O, / . 

const = H-(84) 

r 2 
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With this choice one obtains the following formulas for the self energy which are fully 
symmetric under the transformation N —» —N: 


(V± 2 (v)) = (3^—(v 3 ) 

(vm) = ± 1 -r{-\ {Lp2)+ ^y ) 

Written in terms of the dimensionless held variable $ = p/p m ax this is 

(X4i(^)> = ±^pL±(* 2 ) 


(V^v)) = (|<4> 3 } ~ <*>) 


(V± 3 (^)) = ±- p 2 m 


-(J. 2 ) + (J. 4 ) 


For Tchebyscheff maps of arbitrary order N one obtains 


V±n(<p) = — ~~pLox ± J cos(A^arccos$)d$| 

= Hr*- {-* 2 * (jfTT w * ) - j^T r *-‘ ( * ) ) } 


and 


(85) 

( 86 ) 
(87) 


( 88 ) 

(89) 

(90) 

(91) 

(92) 


(V±*M) = (±1) W F^^ + cj. (93) 

For uncoupled Tchebyscheff maps with \N\ > 2, any expectation of an observable y4(<h) is 
given by eq. m and (plj) . For a ^ 0 the invariant density changes in a nontrivial way, but 
expectations can still be easily calculated numerically by long-time iteration of the coupled 
map lattice. 
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Fig. 3 Self energy V(a) of the type-A chaotic field in the low-coupling region. There are 
local minima at couplings a* that coincide with the weak coupling constants of right-handed 
fermions in the standard model. 




